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PU  RPOSE :  Iron  has  been  used  to  reduce  organic  contaminants  including  high  explosives  (HE)  such 
as  Research  Department  Explosive  (RDX).  Work  is  ongoing  at  the  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC)  to  develop  guidelines  for  using  iron  bed  reactors  to  remove  RDX  in 
surface  water  runoff  from  the  impact  areas  of  firing  ranges.  As  part  of  this  work,  a  mathematical  model 
was  developed  to  gain  a  better  understanding  of  the  experimental  results  and  to  provide  information  for 
design  and  operation  of  these  reactors.  This  model  and  its  application  results  are  described  within  this 
technical  note  (TN). 

Following  additional  validation  against  field  data,  this  model  will  be  implemented  within  the  Training 
Range  Environmental  Evaluation  and  Characterization  System  (TREECS™)  (http://el. erdc. usace. 
armv.mil/treecs/).  TREECS™  was  developed  by  ERDC  for  the  Army  with  varying  levels  of  capability 
to  forecast  the  fate  of  munitions  constituents  (MC),  such  as  HE  and  metals,  within  and  transported  from 
firing/training  ranges  to  surface  water  and  groundwater.  The  overall  purpose  of  TREECS™  is  to 
provide  environmental  specialists  with  tools  to  assess  the  potential  for  MC  migration  into  surface  water 
and  groundwater  systems  and  to  assess  range  management  strategies  to  ensure  protection  of  human 
health  and  the  environment.  Although  TREECS™  was  developed  to  forecast  the  fate  of  MC  on  firing 
ranges,  it  has  general  applicability  to  many  other  situations  requiring  prediction  of  contaminant  fate  in 
multi-media  environmental  systems. 

BACKGROUND:  Munitions  containing  HE  fired  into  impact  areas  can  infrequently  experience  low- 
order  detonations  that  result  in  unexploded  MC  residue.  MC  residues  that  build  up  over  time  can  be 
transported  with  water  to  off-range  receiving  water,  potentially  resulting  in  concentrations  exceeding 
protective  health  guidelines.  Methods  are  sought  that  can  be  implemented  in  the  field  to  remove  RDX 
in  runoff  from  firing  range  impact  areas.  Metallic  iron  can  be  mixed  with  other  materials,  such  as 
gravel,  biodegradable  organic  matter,  and  iron-reducing  bacteria,  to  create  a  reactor  bed  for  reducing 
RDX  and  its  degradation  products.  For  field  conditions,  the  reactor  can  be  placed  within  a  trench  that 
collects  rainfall  runoff  from  the  firing  range  impact  area.  The  reactor  is  water-saturated,  which,  in  the 
presence  of  oxygen,  causes  the  metallic  iron  to  rust  fonning  Fe(II)  and  Fe(III).  Acetate  is  added  to  the 
reactor  providing  food  substrate  for  the  bacteria.  As  the  substrate  is  utilized  by  the  bacteria,  oxygen  and 
Fe(III)  serve  as  electron  acceptors,  thus  reducing  most  of  the  iron  to  Fe(II).  When  RDX  is  introduced 
into  the  reactor,  Fe(II)  serves  as  an  electron  donor,  and  RDX  serves  as  the  electron  acceptor,  thus 
oxidizing  Fe(II)  to  Fe(III)  while  reducing  RDX  through  the  sequence  of  degradation  products  of  MNX, 
DNX,  TNX,  and  eventually  various  end-products,  such  as  formaldehyde.  If  RDX  is  continually 
supplied  to  the  reactor,  much  of  the  Fe(II)  can  be  converted  to  Fe(III),  and  the  reactor  loses  its  ability  to 
degrade  RDX.  Thus,  acetate  must  be  re-supplied  to  reduce  Fe(III)  to  Fe(II).  Additionally,  some  of  the 
Fe(II)  can  become  dissolved  and  flushed  out  by  the  water  flowing  through  the  reactor. 
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A  degradation  pathway  for  RDX  is  shown  in  Figure  1 .  Two  electrons  are  required  for  each  degradation 
step.  The  proposed  reactions  (Szecsody  et  al.  2001)  in  these  degradation  processes  are: 

RDX  +  2Fe2+  <=>  MNX  +  2Fe3+  +  20H' 

MNX  +  2Fe2+  <=>  DNX  +  2Fe3+  +  20H' 

DNX  +  2Fe2+  <=>  TNX  +  2Fe3+  +  20H' 

TNX  +  2Fe2+  <=>  uni  +  2Fe3+  +  20H' 

uni  +  14Fe2+  <=>  CH3NHNHCH3  +  MeOH  +  HCHO  +  14Fe3+  +  140H'  (1) 

where  uni  is  a  non-cyclic  unknown  compound  that  degrades  into  benign  by-products.  Shrout  et  al. 
(2004)  suggest  that  uni  is  three  moles  of  hydroxymethylnitrosamine  following  hydrolysis. 


A  study  by  Wanaratna  et  al.  (2006)  clearly  indicates  that  their  observed  RDX  degradation  rates  were 
dependent  upon  the  amount  of  solid  phase,  zero-valent  iron  (ZVI)  within  their  reactor,  and  their 
results  support  the  theory  that  the  reduction  of  RDX  occurs  on  the  surfaces  of  the  solid  iron.  These 
results  are  important  since  they  indicate  that  the  model  development  should  proceed  along  the  path 
of  surface  site  reaction  following  mass  transfer  of  dissolved  RDX  from  pore- water  to  solid  surfaces. 
This  modeling  approach  is  quite  different  from  an  approach  representing  reactions  of  dissolved  RDX 
and  iron  in  water. 


OBJECTIVE:  The  objective  of  this  study  was  to  develop  a  model  to  predict  the  degradation  of 
RDX  and  its  degradation  products  within  an  iron  bed  reactor.  A  batch  reactor  model  and  a  one¬ 
dimensional  (ID),  longitudinal  reactor  model  are  derived  for  surface  site  reactions  and  tested  against 
laboratory  experiments  as  described  below. 

MODEL  FORMULATION  FOR  BATCH  REACTOR:  Mass  balance  equations  with  first-order 
degradation  kinetics  can  be  used  to  model  the  concentrations  of  RDX  and  its  degradation  products 
shown  in  Equation  set  1 .  These  equations  are  first  developed  for  a  single  fully  mixed  (batch)  reactor 
consisting  of  a  water-saturated  porous  media,  with  a  constant  RDX  loading,  and  for  a  constant  water 
flow  rate  through  the  reactor.  These  equations  are  stated  as  follows: 

dRDX  =  RDX^Q_RDXQ_  x 
dt  <j)VT  <t>VT 


2 


d-M^  =  KRDX-^Q 

dt  (j)VT 

d  DNX  ,  t/fATV  DNXQ 
dt  (f)VT 

dJ^L=KDNX.mQ 


iJ^.=kjNx-u^Q 

dt  (j)VT 


-  k2  MNX 

-  k3  DNX 

-  k4  TNX 

k5  uni 
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(3) 

(4) 

(5) 

(6) 


where 


ki, 2, 3. 4.5 


RDXin 

Q 

t 

<t> 

VT 


overall  degradation  rates  including  mass  transfer  to  surface  sites  and  surface  site 
reaction  rate,  hr’1 

concentration  of  RDX  flowing  into  the  reactor,  mg/L 

flow  rate  of  water  through  the  reactor  and  carrying  RDX  into  the  reactor,  L/hr 
time,  hr 

porosity  of  the  reactor  media,  water  volume/total  volume 
reactor  total  volume,  liters  (L) 


Concentrations  of  all  constituents  (RDX,  MNX,  etc.)  are  expressed  in  mg/L.  The  product  RDXinQ  is  the 
RDX  loading  rate  (mg/hr)  into  the  reactor.  It  is  assumed  that  all  constituents  (RDX  and  degradation 
products)  are  dissolved  in  water  (i.e.,  mass/volume  water)  and  that  sorption  partitioning  is  minor  and 
can  be  ignored.  This  last  assumption  is  reasonable  given  the  range  of  partitioning  coefficients  for  RDX 
and  its  degradation  products.  It  is  assumed  for  this  model  that  the  concentration  of  Fe(II)  is  plentiful 
and  continuous  over  time,  which  means  that  acetate  or  a  similar  food  source  for  reducing  bacteria  is 
regularly  provided  to  the  reactor.  Therefore,  this  model  does  not  attempt  to  compute  the  supply  or 
concentration  of  Fe(II). 

For  this  model  development,  it  is  assumed  that  all  of  the  RDX  degradation  products  have  the  same 
reaction  rate  as  RDX,  i.e.,  kj=k2=k3=k4=k5=k.  As  mentioned  above,  the  reaction  rate  k  (hr'1)  is  an 
overall  reaction  rate  that  accounts  for  mass  transfer  from  pore  water  to  solid  surfaces  and  surface-site 
reaction  (i.e.,  reduction).  The  immediate  goal  is  to  derive  a  relationship  for  detennining  k. 

Mass  transfer  here  refers  to  the  movement  of  solute  (e.g.,  RDX)  from  pore  water  to  solid  surface 
(iron).  For  a  general  solute  concentration  C  (mg/L),  this  mass  transfer  can  be  described  as  follows: 

^  =  -Ka,(C-C„)  (7) 

where 


km  =  the  mass  transfer  coefficient,  m/hr'1 
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as  =  specific  surface  area  of  solids,  which  is  the  solid  surface  area  a  divided  by  the  pore- 
water  volume  Vw,  m'1 

Cb  =  solute  concentration  at  the  solid  surface,  mg/L 

The  problem  with  Equation  7  is  that  the  parameters  km,  as,  and  G>  are  difficult  to  determine.  To 
simplify  matters,  the  product  of  the  first  two  parameters  can  be  represented  by  the  mass  transfer  rate 
X  (per  hour),  or 


dC_ 

dt 


yc-c„) 


(8) 


An  approach  for  estimating  X  is  presented  later  below. 

The  solute  concentration  at  the  surface  of  the  solid  Q,  is  driven  by  the  surface-site  reaction 


dt 


KC, 


(9) 


where  K  is  the  surface-site  reaction  rate  (per  time).  Wanaratna  et  al.  (2006)  tested  RDX  reduction 
using  zero-valent  iron  (ZVI)  and  detennined  that  K  is  related  to  the  amount  of  ZVI  through  a  Monod 
relationship 


K  Kjzvi] 

K+\zvi] 


(10) 


with  kmax  =  2.824  min"1  or  0.047  sec"1,  and  ko  =  1.806  M  (mole/L).  The  concentration  of  ZVI  is  also  in 
moles/L  and  can  be  detennined  by  dividing  the  total  mass  (g)  of  ZVI  in  the  reactor  by  the  reactor  water 
volume  {(f>VT)  and  dividing  the  result  by  the  molecular  weight  of  iron  (55.85  g/mole).  For  the  present 

reactor  conditions,  it  is  the  concentration  [Fe(II)]  that  is  required  in  Equation  10  rather  than  [ZVI]. 
Thus,  the  values  for  kmax  and  k0  in  Equation  10  could  be  different  for  Fe(II)  than  for  ZVI.  If  the 
concentration  of  Fe(II)  is  much  greater  than  the  half-saturation  concentration  for  Fe(II),  i.e.,  ko,  then  K 
is  simply  equal  to  Gar- 


Considering  that  equilibrium  conditions  should  eventually  exist  between  mass  transfer  and  surface- 
site  reaction,  Equations  8  and  9  can  be  set  equal  to  each  other,  resulting  in 


X(C-Cb)  =  KC„ 


(11) 


or 


Q 


xc 

X  +  K 


Substituting  the  above  relation  for  Q,  into  Equation  8  yields 


(12) 


4 


dC_ 

dt 


XK 


X  +  K 


C 


-kC 
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(13) 


Thus,  the  overall  reaction  rate  k  accounting  for  mass  transfer  and  surface-site  reaction  is 


,  XK 

k  — - 

X  +  K 


(14) 


When  K  is  large  compared  to  A,  the  overall  reaction  rate  is  driven  by  mass  transfer.  When  A  is  large 
compared  to  K,  the  overall  reaction  rate  is  driven  by  the  surface-site  reaction  rate. 

Mass  transfer  between  fluid  and  liquid  and  fluid  and  solid  has  been  extensively  studied  over  the  past 
50  years.  Many  such  studies  involved  dissolving  spheres  or  cylinders  in  packed  beds  or  suspended  in 
fluids.  Although  various  correlations  have  been  used  to  describe  theory  and  experimental  results, 
there  is  no  one  single  correlation  that  seems  to  work  for  all  conditions.  Most  of  the  correlations 
relate  the  Sherwood  number  (Sh)  to  the  Peclet  number  (Pe). 

Many  such  correlations  were  reviewed  and  tested  against  the  meso-scale  laboratory  results  described 
below  in  this  TN.  The  one  correlation  that  appeared  to  hold  the  most  promise  is  based  on  data 
presented  by  Miyauchi  (1971).  Figure  6  within  the  paper  by  Miyauchi  (1971)  shows  a  strong 
relationship  between  Sh  and  Pe,  where  Sh  is  defined  as  kmdp/De  and  km  is  the  mass  transfer  coefficient 
(m/sec);  dp  is  the  diameter  (m)  of  the  packed  bed  material;  and  De  is  the  effective  diffusivity  (m2/sec)  of 
the  solute  in  fluid  within  the  pore  spaces.  The  Peclet  number  is  defined  as  udp/De,  where  u  is  the  pore 
velocity  (m/sec).  The  pore  velocity  is  the  Darcy  velocity  U divided  by  the  bed  porosity,  and  the  Darcy 
velocity  is  the  reactor  water  flow  rate  divided  by  the  cross-sectional  area  of  the  flow.  For  practical 
purposes,  it  is  assumed  that  De  is  equivalent  to  the  molecular  diffusivity  Dm.  An  approximate  fit  of  the 
data  presented  in  Figure  6  of  Miyauchi  (1971)  resulted  in  the  relationship 

Sh  =  APem  (15) 


Similar  relationships  have  been  reported  by  others,  including  one  documented  by  Wilson  and 
Geankoplis  (1966),  which  was  derived  from  the  dissolution  of  solid  benzoic  acid  spheres. This 
relationship  was  also  used  successfully  by  Geller  and  Hunt  (1993)  to  predict  the  mass  transfer  from 
non-aqueous  phase  organic  liquids  to  water  in  porous  media.  However,  the  relationship  described  by 
Wilson  and  Geankoplis  (1966)  was  not  developed  for  dilute  beds  and  under-predicted  the  meso-scale 
results.  Miyauchi  (1971)  presented  data  for  dilute  beds,  where  the  reactive  particles  were  rather  sparse 
relative  to  the  packing  particles.  This  is  the  situation  for  solid  iron  placed  in  gravel  beds,  which  is  the 
case  for  the  meso-scale  experiments  and  will  most  likely  be  the  case  for  field-scale  reactors. 

Substituting  definitions  for  Sh  and  Pe,  Equation  15  can  be  solved  for  km 


1/3  r-v  2/3 


k,„  = 


4  uUiD. 


d 


2/3 


(16) 
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Assuming  that  the  packing  particles  are  spherical,  a  relationship  for  specific  surface  area  as  can  be 
derived  from  the  relationships  for  the  surface  area  ( a )  of  an  individual  particle  and  volume  of  water 
(Vw)  along  with  the  definition  of  void  fraction  or  porosity.  This  relationship  is 


a  _  6(1-^) 
V  d  d) 

W  p  T 


(17) 


where  (f>  is  the  porosity  as  before.  Similarly,  a  relationship  for  cylindrical  packing  particles  can  be 
derived 


a 


s 


2(1  —  ^) 


(18) 


where  dp  is  the  diameter  (m)  of  an  individual  cylindrical  particle,  and  /  is  its  length  (m).  However,  in 
most  cases,  it  can  be  assumed  that  the  packing  particles  are  spherical.  The  mass  transfer  rate  X  (sec'1) 
for  spherical  packing  particles  can  now  be  obtained 


X 


24  umD2'3 


(w) 


5/3 


(19) 


For  cylindrical  packing  particles,  the  mass  transfer  rate  is 


2 


k  a  — 

m  s 


8  «i1/3Af(W) 


HP 


5/3 


d 

— +  2 
/ 


(20) 


When  the  cylinder  length  is  equal  to  the  cylinder  diameter,  Equation  20  becomes  the  same  as 
Equation  19. 

It  is  necessary  to  allow  for  a  minimum  mass  transfer  rate  for  no  flow  or  stagnant  pore-water 
conditions.  Even  under  stagnant  conditions,  there  will  be  some  solute  movement  due  to  molecular 
diffusion.  The  literature,  including  Miyauchi  (1971),  consistently  indicated  that  a  minimum  Sh  of 
approximately  2.0  describes  mass  transfer  when  Pe  «  1 ,  or  when  there  is  essentially  no  flow.  Thus, 
Sh  =  2  was  used  for  setting  a  minimum  mass  transfer  rate,  resulting  in  the  following  relationship  for 
spherical  packing  particles: 


no.  (i-<*) 
K 


The  minimum  mass  transfer  rate  for  cylindrical  packing  particles  is 


4,in  - 


4  A„  (!-«>) 

\dp 

1 

\ 

(21) 


(22) 
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Using  either  Equation  19  or  20  (with  Equations  21  and  22  for  lower  limits)  and  an  estimate  of  K 
from  Equation  10,  A:  can  be  obtained  from  Equation  14,  and  Equations  2-6  can  be  solved  for  RDX 
and  its  degradation  products.  It  is  assumed  that  at  time  zero,  the  RDX  concentration  in  the  reactor  is 
known,  and  the  initial  concentrations  of  MNX,  DNX,  TNX,  and  uni  are  zero.  The  flow  rate  and 
RDX  concentration  of  water  entering  the  reactor  are  also  known,  as  well  as  the  mass  of  iron  in  the 
reactor,  which  is  assumed  to  be  Fe(II).  With  the  input  parameters,  initial  conditions,  and  boundary 
conditions,  the  five  differential  equations  are  solved  in  a  spreadsheet  yielding  the  constituent 
concentration  over  time,  where  the  time-step  of  the  solution  is  specified  by  the  user.  The  first-order- 
accurate,  Euler  integration  method  is  used  for  the  solution. 

A  conservative  tracer  is  also  included  within  the  model  to  evaluate  flushing  time.  The  mass  balance 
equation  for  a  tracer  moving  through  the  fully  mixed  reactor  is 


dC 

dt 


Q_ 

<j>VT 


K-c) 


(23) 


where  C  is  the  tracer  concentration  in  the  reactor,  and  C„,  is  the  tracer  concentration  in  the  inflow. 

SOLUTIONS  WITH  BATCH  REACTOR  MODEL:  The  batch  reactor  model  was  first  applied  to  a 
hypothetical  reactor  for  several  conditions  to  gain  an  understanding  of  basic  model  behavior.  The 
first  test  case  consisted  of  no  water  flow  and  no  RDX  loading  into  the  reactor  but  with  an  initial 
RDX  concentration  within  the  reactor.  Reactor  input  conditions  are  shown  in  Table  1.  The  two 
reaction  kinetic  parameters  were  set  to  those  detennined  for  ZVI  by  Wanaratna  et  al.  (2006)  pending 
availability  of  additional  data  for  Fe(II)  as  the  reactant. 


Table  1.  Input  parameters  for  hypothetical  batch  reactor  model,  first  test  case. 

Input  Parameter 

Input  Value 

Units 

Reactor  volume,  VT 

100 

L 

Media  porosity,  <f> 

0.5 

dimensionless 

Water  flow  rate  entering/exiting,  Q 

0 

L/hr 

RDX  concentration  flowing  in,  RDXin 

0 

mg/L 

Initial  RDX  concentration  in  reactor 

1.0 

mg/L 

Mass  of  iron  in  the  reactor 

1000 

g 

Average  particle  diameter  of  packed  bed  material,  dp 

.01 

m 

Molecular  diffusivity  of  RDX  and  degradation  products,  Dm 

5.9E-10 

m2/sec 

Maximum  surface  site  reaction  rate,  kmax 

2.824 

min"1 

Half-saturation  concentration  of  iron  for  surface  site  reaction,  k0 

1.806 

mole/L 

The  time  series  results  for  this  test  case  are  shown  in  Figure  2  for  RDX  and  its  degradation  products. 
Since  there  is  no  loading  of  RDX  after  the  initial  specified  concentration,  the  RDX  concentration 
drops  over  time.  The  degradation  products  build  up  and  then  die  off  over  time  in  the  sequence  of  the 
pathway  shown  in  Figure  1. 
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Figure  2.  Concentrations  of  RDX  and  degradation  products 
versus  time  for  hypothetical  batch  reactor,  first  test 
case. 


The  second  test  for  the  hypothetical  reactor  was  for  a  constant  water  flow  rate  of  10  L/hr,  but  with 
zero  RDX  in  the  inflow.  All  other  inputs  remained  the  same  as  the  first  test  case,  which  included  an 
initial  RDX  concentration  of  1.0  mg/L.  The  inputs  for  this  test  condition  are  shown  in  Table  2. 


Table  2.  Input  parameters  for  hypothetical  batch  reactor  model,  second  test  case. 

Input  Parameter 

Input  Value 

Units 

Reactor  volume,  VT 

100 

L 

Media  porosity, 

0.5 

dimensionless 

Water  flow  rate  entering/exiting,  Q 

10 

L/hr 

RDX  concentration  flowing  in,  RDXin 

0 

mg/L 

Initial  RDX  concentration  in  reactor 

1.0 

mg/L 

Mass  of  iron  in  the  reactor 

1000 

g 

Average  particle  diameter  of  packed  bed  material,  dp 

.01 

m 

Molecular  diffusivity  of  RDX  and  degradation  products,  Dm 

5.9E-10 

m2/sec 

Maximum  surface  site  reaction  rate,  kmax 

2.824 

min"1 

Half-saturation  concentration  of  iron  for  surface  site  reaction,  k0 

1.806 

mole/L 

The  results  of  the  second  test  are  shown  in  Figure  3  for  RDX  and  its  degradation  products.  These 
results  are  similar  to  those  of  the  first  test  case,  except  that  the  decline  of  concentrations  is  faster  due  to 
flushing  with  clean  water.  A  tracer  concentration  of  100  mg/L  was  introduced  within  the  reactor 
influent.  The  tracer  concentration  versus  time  for  the  second  test  is  shown  in  Figure  4.  The  tracer 
approaches  the  inflowing  tracer  concentration  of  1 00  mg/L  as  time  progresses.  The  hydraulic  residence 
time  of  the  reactor  is  5  hr;  thus,  it  takes  longer  than  this  for  the  reactor  to  reach  the  steady-state 
concentration  of  100  mg/L.  The  tracer  results  agree  with  analytical  results  for  this  test  case,  such  as 
requiring  1 1 .5  hr  to  reach  a  level  of  90%  completely  flushed. 

The  third  test  case  for  the  hypothetical  batch  reactor  was  conducted  with  a  constant  RDX  loading  as 
prescribed  for  a  water  flow  rate  of  10  L/hr  with  an  RDX  concentration  of  1.0  mg/L  in  this  flow.  The 
initial  RDX  concentration  in  the  reactor  was  set  to  zero.  All  other  inputs  were  the  same  as  for  the  first 
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Figure  3.  Concentrations  of  RDX  and  degradation  products 
versus  time  for  hypothetical  batch  reactor, 
second  test  case. 


120 


Time,  hrs 


Figure  4.  Tracer  concentration  versus  time  for  hypothetical 
batch  reactor,  second  test  case. 


Table  3.  Input  parameters  for  hypothetical  batch  reactor  model,  third  test  case. 

Input  Parameter 

Input  Value 

Units 

Reactor  volume,  VT 

100 

L 

Media  porosity,  </> 

0.5 

dimensionless 

Water  flow  rate  entering/exiting,  Q 

10 

L/hr 

RDX  concentration  flowing  in,  RDXin 

1.0 

mg/L 

Initial  RDX  concentration  in  reactor 

0 

mg/L 

Mass  of  iron  in  the  reactor 

1000 

g 

Average  particle  diameter  of  packed  bed  material,  dp 

.01 

m 

Molecular  diffusivity  of  RDX  and  degradation  products,  Dm 

5.9E-10 

m2/sec 

Maximum  surface  site  reaction  rate,  kmax 

2.824 

min"1 

Half-saturation  concentration  of  iron  for  surface  site  reaction,  k0 

1.806 

mole/L 
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two  tests.  The  inputs  for  this  test  case  are  shown  in  Table  3.  The  results  of  this  test  case  are  shown  in 
Figure  5 .  RDX  reaches  a  steady-state  concentration  balanced  by  losses  due  to  degradation  and  flushing 
and  gains  due  to  the  constant  loading  of  RDX.  Likewise,  constant  concentrations  of  RDX  degradation 
products  are  produced.  The  results  for  all  three  of  the  hypothetical  batch  reactor  tests  appear  to  be 
reasonable  and  expected. 


-  RDX 
■  MNX 
*  DNX 
TNX 
uni 


Figure  5.  Concentrations  of  RDX  and  degradation 

products  versus  time  for  hypothetical  batch 
reactor,  third  test  case. 

The  batch  reactor  model  was  next  applied  to  laboratory  tests  conducted  in  the  ERDC  Environmental 
Laboratory.  These  tests  were  conducted  with  100-mL  vials  containing  iron  filings  and  spiked  with 
soil  as  inoculants  and  an  initial  dissolved  RDX  concentration  of  4.0  mg/L.  A  metallic  iron  mass  of 
0.02  g  was  added  to  the  vials  initially.  Thus,  each  vial  contained  3.58  mmol/L  of  iron.  The  RDX 
concentration  in  the  vials  was  measured  over  60  days.  Tests  were  conducted  at  30  and  13.5  °C  in 
MaxQ  6000  Incubated/Refrigerated  Shakers  (Thermo  Scientific,  USA).  The  measured  RDX 
concentrations  decreased  with  increasing  time,  but  there  was  little  difference  in  results  for  the  two 
temperature  conditions.  Therefore,  the  observed  concentrations  for  the  two  temperatures  were 
averaged  for  each  sampling  time  point.  The  vials  were  shaken  continuously  at  130  RPM  during  the 
experiment.  The  iron  filings  were  approximately  0.5  mm  wide  and  1  mm  long.  Thus,  the  filings  were 
assumed  to  be  cylindrical. 

The  porosity  was  computed  from 


^1-^1-M 


Lr 


PsVt 


(24) 


3  3 

where  Vs  is  the  volume  (m  )  of  solid  iron,  Vt  is  the  reactor  volume  (m  ),  Ms  is  the  mass  (g)  of  solid 
iron  in  the  system,  and  ps  is  the  density  of  solid  iron  of  7.87E6  g/m  .  It  is  not  known  how  fast  the 
water  was  moving  within  the  vials  when  shaken,  so  u  (pore  velocity)  was  varied  until  model  results 
matched  observed  data.  Using  a  value  of  zero  for  u  resulted  in  insufficient  RDX  reduction  over  time. 
The  calibrated  value  of  u  was  1.0E-5  m/sec.  The  values  for  the  two  kinetic  parameters  were  set  to 
the  values  used  for  the  previous  tests.  All  input  values  for  this  test  case  are  summarized  in  Table  4. 
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Table  4.  Input  parameters  for  batch  reactor  model,  laboratory  experiment  with  100-mL 
vial. 

Input  Parameter 

Input  Value 

Units 

Reactor  volume,  VT 

0.1 

L 

Media  porosity,  (j> 

0.999975 

dimensionless 

Water  flow  rate  entering/exiting,  Q 

0 

L/hr 

RDX  concentration  flowing  in,  RDXin 

0 

mg/L 

Initial  RDX  concentration  in  reactor 

4 

mg/L 

Mass  of  iron  in  the  reactor 

0.02 

g 

Average  particle  diameter  of  packed  bed  material,  dp 

.0005 

m 

Average  particle  length  of  packed  bed  material,  / 

.001 

m 

Molecular  diffusivity  of  RDX  and  degradation  products,  Dm 

5.9E-10 

m2/sec 

Maximum  surface  site  reaction  rate,  kmax 

2.824 

min"1 

Half-saturation  concentration  of  iron  for  surface  site  reaction,  k0 

1.806 

mole/L 

Calibrated  u 

1.0E-5 

m/sec 

Observed  and  model-computed  RDX  concentrations  for  the  vial  experiment  are  compared  in  Figure  6. 
The  model  matches  the  observed  data  closely,  but  it  was  necessary  to  adjust  u  to  reach  such  agreement 
since  the  value  of  this  variable  was  not  known  for  shaken  vials.  A  greater  value  of  u  causes  under¬ 
prediction  of  concentration,  and  a  smaller  value  causes  over-prediction.  This  application  was  not  an 
adequate  validation  of  the  model  due  to  the  need  to  adjust  u. 


Figure  6.  RDX  concentration  versus  time  for  laboratory  vial 
experiment. 

MODEL  FORMULATION  FOR  ID  LONGITUDINAL  REACTOR:  This  model  is  designed  to 
account  for  constituent  variations  along  the  longitudinal  axis  of  the  reactor  as  well  as  transport 
through  the  mobile  domain  and  constituent  storage  in  the  immobile  domain.  Consider  a  water- 
saturated,  longitudinal  reactor  of  uniform  width,  depth,  and  porous  media  properties  (e.g.,  porosity 
(/>)  that  receives  a  water  flow  rate  and  loading  of  RDX  at  the  upstream  end.  For  a  given  water  inflow 
rate  Q  (L/hr),  reactor  width  W (m),  and  depth  H(m),  the  Darcy  velocity  U  (m/hr)  through  the  reactor 

can  be  computed  from  U  = - — - .  Dispersion  can  also  be  included  in  the  transport  along  the 

1000  WH 
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reactor.  Transport  trapping  is  included  through  the  use  of  a  dual-domain  modeling  approach  (Zheng 
and  Wang  1999,  Liu  et  al.  2007).  Dual-domain  models  provide  a  means  to  represent  the  effects  of 
dead-end  pores.  The  time-varying,  one-dimensional  reactive  transport  equation  for  the  longitudinal 
reactor  for  mobile  and  immobile  solute  constituent  concentrations  is 


—  (C.  -C  )  +  D 

r  i  \  iyn  m )  1  a: 


=  ~  ~  d2c.„  U  dc,„ 

dt  ftp 

dc, 


dx  ftp  dx 


- kC+S , 


m  '  m 


dt  (1  -f)<p 


[c  -C  )-k  .  c  +  s. 

\  m  im )  min  im  1  in 


(25) 

(26) 


where 


Cm 

Cim 


f  = 


X  = 

t  = 
Dx  = 
k  = 

' min 


u  = 


solute  concentration  in  the  mobile  domain  (mg/L) 

solute  concentration  in  the  immobile  domain  (mg/L) 

fraction  of  domain  that  is  mobile  or  tpm/tp  where  <pm  is  the  mobile  or  effective 

porosity 

distance  along  the  longitudinal  axis  of  the  reactor  (m) 
time  (hr) 

dispersion  coefficient  (mVhr) 

overall  reaction  rate  in  the  mobile  domain  accounting  for  mass  transfer  and  surface- 
site  reaction  (1/hr) 

the  minimum  overall  reaction  rate  accounting  for  mass  transfer  and  surface-site 


reaction  for  no  flow  in  the  immobile  domain  (1/hr),  kmm 


K  mK 


+  K 


Darcy  velocity  (m/hr) 

bed  porosity  (dimensionless) 

source  term  in  the  mobile  domain  (mg/L/hr) 

source  term  in  the  immobile  domain  (mg/L/hr) 

first-order  mass  transfer  rate  between  mobile  and  immobile  domains  (1/hr) 


The  variable  k  is  the  same  as  previously  defined  for  the  batch  reactor.  It  is  noted  that  the  two  above 
equations  do  not  include  sorption  partitioning;  thus,  the  solute  is  assumed  to  have  a  very  low 
partitioning  coefficient,  which  is  generally  the  case  for  RDX  and  its  degradation  products. 

The  source  terms  Sm  and  Sun  are  zero  for  RDX.  The  source  tenn  for  MNX  is  the  kinetic  loss  tenn  for 
RDX,  or 


Sm(MNXm) 


XK 
X  +  K 


RDX 

m 


(27) 
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S„(MNX,„)  =  ^-RDX,m  (28) 

Likewise,  the  source  tenn  for  DNX  is  the  kinetic  loss  term  for  MNX,  and  so  on  down  the  reaction 
chain.  This  model  only  allows  for  spherical  bed  solids.  The  mass  transfer  rate  for  the  mobile  domain 
is  set  to  the  maximum  value  of  A  computed  from  Equation  1 9  and  Amm  computed  from  Equation  2 1 , 
thus  allowing  for  cases  when  the  flow  velocity  is  small.  The  mass  transfer  rate  for  the  immobile 
domain  is  Amm  as  determined  from  Equation  2 1 . 

Dispersivity  is  used  within  the  model  to  compute  the  spatially  varying,  longitudinal  dispersion 
coefficient  as  follows: 


Dx  =ax  XU  (29) 

where  axX is  the  longitudinal  dispersivity,  ax  is  the  proportionality  constant  (dimensionless),  andX 
is  the  distance  from  the  beginning  of  the  reactor  (where  water  flows  in)  to  the  location  of  interest.  A 
typical,  commonly  used  value  of  the  longitudinal  dispersivity  constant  is  0.1. 

Equations  25  and  26  include  two  unknowns  (Cm  and  Cim)  that  must  be  numerically  solved  over  the 
two  independent  variables  x  and  t,  the  space  and  time  domains.  The  initial  conditions  for  this  model 
are  zero  concentration  for  all  constituents.  The  boundary  conditions  that  are  specified  atx  =  0  are 
zero  concentration  for  all  constituents  except  for  RDX,  which  is  specified.  The  flow  rate  of  water 
entering  the  reactor,  and  thus  the  flow  speed,  as  well  as  the  concentration  of  RDX  in  the  influent 
water,  can  be  varied  over  time  for  this  more  general-purpose  model.  However,  this  added  flexibility 
requires  a  numerical  solution  that  is  better  suited  to  programming  with  a  computer  coding  language. 

The  above  fonnulations  were  programmed  into  a  numerical  model  code  using  Fortran.  The  finite 
difference  method  was  used  to  discretize  time  and  space  and  solve  the  two  partial  differential 
equations.  The  equations  are  solved  explicitly  using  a  forward  difference  in  time  and  options  of  either 
central  or  upwind  differencing  for  advection  (the  third  tenn  on  the  right-hand  side  of  Equation  25).  A 
central  difference  was  used  for  the  dispersion  tenn  (the  second  tenn  on  the  right-hand  side  of  Equation 
25).  The  model  code  checks  for  numerical  stability  requirements  associated  with  the  solution  schemes. 
It  is  noted  that  the  dispersion  tenn  cannot  be  zero  if  stability  is  to  be  maintained  when  using  the  central 
difference  for  advection. 

SOLUTIONS  WITH  THE  ID  LONGITUDINAL  REACTOR  MODEL:  Hypothetical  reactor 
conditions  closely  similar  to  the  third  test  conditions  described  in  the  batch  reactor  model  applications 
were  applied  to  the  longitudinal  reactor  model.  These  test  conditions  are  shown  in  Figure  7.  The  length 
of  the  reactor  was  set  to  1.0  m,  and  the  width  and  height  were  set  to  0.32  m,  which  yielded 
approximately  the  same  total  volume  of  100  L  as  used  for  the  batch  reactor  test  case.  The  fraction  of 
mobile  domain  was  set  to  1 .0. 

Most  of  the  inputs  shown  in  Figure  7  are  self-explanatory,  but  a  few  require  further  definition.  The 
input  variable  labeled  “domain  mass  trans  rate”  is  ((and  is  shown  in  Equations  25  and  26.  The  variable 
“number  of  delta  x”  is  the  number  of  computational  segments  NX  that  the  total  reactor  length  is  to  be 
discretized.  For  the  example  shown,  AX=10;  thus,  the  spatial  size  of  each  computational  segment  is 
0.1  m.  The  variable  “DX  output  intervals”  prescribes  the  number  of  segments  to  skip  within  the 
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discretized  X  space  for  writing  output.  For  example,  a  value  of  2  means  that  output  will  be  written  for 
every  second  segment,  starting  with  segment  1  and  ending  at  the  last  segment.  The  variable  “solution 
scheme”  should  have  a  value  of  either  0  or  1 .  A  value  of  0  denotes  that  the  upwind  differencing  scheme 
is  used  for  advection,  whereas  a  value  of  1  denotes  that  the  central  differencing  scheme  is  used. 

As  with  the  batch  reactor  model,  two  iron  kinetic  reaction  parameters  must  be  specified,  potentially 
requiring  calibration  against  laboratory  or  field  results.  The  default  values  are  shown  in  Figure  7.  Two 
additional  parameters  (f  and  Q  are  introduced  by  using  the  dual-domain  modeling  approach  to  account 
for  transport  trapping.  Additionally,  the  longitudinal  dispersivity  constant  must  be  specified,  but  the 
default  value  shown  in  Figure  7  is  typically  used.  The  time-step  and  number  of  spatial  segments  may 
have  to  be  adjusted  to  provide  numerical  stability  and  sufficient  accuracy. 


"Test  case  for  generic  ID  reactor  model" 

"Length,  m" 

1.0 

"Width,  m" 

.32 

"Height,  m" 

.32 

"Porosity" 

0.50 

"fraction  mobile" 

1.0 

"Dispersivity  constant" 

0.1 

"domain  mass  trans  rate,  1/hr" 

1.0 

"RDX  mole,  diffus.,  cm2/sec" 

5.9E-6 

"packing  solids  diam,  m" 

.01 

"Fe  mass  in  reactor,  kg" 

1. 

"Fe  half  sat  cone,  kO,  mol/L" 

1.806 

"Fe  max  reac  rate,  kmax,  l/min" 

2.824 

"RDX  initial,  mg/L" 

0.0 

"number  of  delta  x" 

10 

"delta  t,  hr" 

.01 

"simulation  duration,  hr" 

20. 

"DX  output  intervals" 

2 

"solution  scheme" 

1 

"Time,  hr  flowin,  L/hr  RDX  in. 

mg/L  Tracer  in,  mg/L" 

0.0  10.0 

1.0 

1.0 

2.0  10.0 

1.0 

1.0 

Figure  7.  Input  parameters  for  longitudinal  reactor 
model,  hypothetical  test  conditions. 

The  pore  water  concentrations  of  RDX  and  its  degradation  products  in  the  mobile  domain  (as 
computed  by  the  longitudinal  reactor  model)  are  plotted  versus  time  in  Figure  8  for  a  location  0.25  m 
down-gradient  of  the  influent.  All  concentrations  reach  steadystate  after  a  few  hours.  At  this  location, 
the  reduction  of  RDX  has  resulted  in  a  buildup  of  degradation  products  that  exceed  the  RDX 
concentration,  which  is  an  inverse  relationship  compared  to  the  results  shown  in  Figure  5.  Thus,  a 
longitudinal  reactor  can  respond  quite  differently  from  a  batch  reactor. 

The  concentrations  versus  distance  values  for  this  test  condition  after  20  hr  are  plotted  in  Figure  9.  This 
figure  shows  the  rapid  loss  of  RDX  followed  by  the  growth  of  RDX  products  and  their  subsequent 
degradation  with  distance  down-gradient.  The  cross-over  point  at  which  the  RDX  concentration  falls 
below  the  product  concentrations  occurs  at  approximately  X  =  0.2  m. 

The  hypothetical  longitudinal  reactor  test  case  described  above  was  repeated;  however,  a  conservative 
tracer  pulse  (as  shown  in  Figure  10)  was  introduced  in  the  influent  to  evaluate  the  effects  of  including 
an  immobile  domain.  Additionally,  the  spatial  resolution  was  doubled  using  20  segments  rather  than 


14 


ERDC  TN-EQT-13-1 
September  2013 


10.  Otherwise,  the  input  conditions  of  Figure  7  were  imposed.  Three  tests  were  run.  The  first  test  was 
for  a  100%  mobile  domain.  The  second  and  third  tests  were  run  with  a  50%  mobile  domain  for  both 
tests,  but  with  domain  mass  transfer  rates  of  1 .0/hr  and  0. 1/hr.  The  tracer  concentration  results  for  the 
square  tracer  input  pulse  are  shown  in  Figure  1 1 .  The  results  for  100%  mobile  domain  demonstrate  that 
roughly  a  square  distribution  of  tracer  concentration  versus  time  occurs  at  the  sampling  location  (X  = 
0.425  m).  There  is  some  longitudinal  dispersion  in  the  solution  due  to  the  value  set  for  dispersivity  as 
well  as  due  to  some  numerical  diffusion  associated  with  the  solution  scheme.  The  results  for  a  50% 
mobile  domain  and  a  mass  transfer  rate  of  1.0/hr  demonstrate  a  time-delayed  and  more  diffuse  tracer 
distribution.  Higher  domain  mass  transfer  rates  showed  little  change  from  this  result  for  these  test 
conditions.  The  results  for  a  50%  mobile  domain  and  a  mass  transfer  rate  of  0. 1/hr  follow  those  of  the 
100%  mobile  domain  initially,  but  the  peak  concentration  is  reached  later,  and  the  tailing 
concentrations  recede  more  gradually.  Overall,  these  test  results  appear  to  follow  logical 
considerations. 
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Figure  8.  Concentrations  of  RDX  and  degradation 
products  versus  time  at  X  =  0.25  m  for 
hypothetical,  longitudinal  reactor  test  case. 


Figure  9.  Concentrations  of  RDX  and  degradation 

products  versus  X  after  20  hr  for  hypothetical, 
longitudinal  reactor  test  case. 
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Time,  hrs 


Figure  10. Input  tracer  pulse  entering  hypothetical 
longitudinal  reactor. 


X  =  0.425  m 


♦  100%  mobile 


■  50%  Mobile,  Rate  =  1.0  per 
hr 

A  50%  Mobile,  Rate  =  0.1  per 
hr 


Figure  1 1 .  T racer  concentration  versus  time  for  tracer 
pulse  in  longitudinal  reactor  at  location  X  = 

0.425  m  with  100%  and  50%  mobile  domain. 

The  ID,  longitudinal  reactor  model  was  applied  to  meso-scale  laboratory  tests  conducted  in  the  ERDC 
Enviromnental  Laboratory.  The  meso-scale  reactor  consisted  of  a  4-ft-long  by  1  -ft,  4-in.-high  by  1  -ft, 
8-in.-wide  chamber  filled  with  a  porous  media  consisting  of  gravel,  cypress  mulch,  and  iron  filings. 
The  porous  media  was  saturated  with  water.  Acetate  was  added  to  the  system  to  specifically  feed  the 
iron-reducing  organisms  (Madigan  et  al.  2002)  naturally  occurring  on  the  reactor  packing.  After  the 
iron-reducing  microbes  established  a  reducing  enviromnent,  dissolved  RDX  was  introduced  in  the 
water  influent,  which  had  a  constant  flow  rate  of  about  120  L/hr.  The  RDX  concentration  in  the 
influent  was  relatively  constant  at  about  0.77  mg/L.  The  chamber  was  baffled  into  four  compartments 
where  the  flow  moved  vertically  downward  within  the  first  compartment  passing  underneath  the  first 
baffle,  and  then  the  flow  moved  vertically  upward  in  the  second  compartment  before  passing  over  the 
second  baffle.  Flow  passed  under  the  third  baffle  and  then  over  the  end  wall  as  it  exited  the  chamber.  A 
sampling  port  was  placed  in  each  chamber  compartment  with  the  intake  located  8  in.  below  the  water 
surface.  Thus,  the  first  sampling  port  was  8  in.  (0.2  m)  from  where  the  influent  entered  the  chamber. 
Pore- water  concentrations  were  measured  at  each  sampling  port. 
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Multiple  laboratory  tests  were  conducted  with  the  meso-scale  reactor  using  approximately  the  same 
test  conditions.  Multiple  tests  were  run  since  it  was  difficult  to  maintain  a  constant  flow  and  input 
RDX  concentration.  The  results  from  tests  2,  3,  and  4  were  averaged  to  obtain  a  single  influent  RDX 
concentration  and  time-varying,  observed  RDX  concentrations  at  the  first  sampling  port.  RDX  was 
not  detected  at  the  other  sampling  ports.  RDX  degradation  products  were  not  measured.  With  a 
constant  influent  RDX  loading  rate  and  water  flow  rate,  steady-state  conditions  should  be  reached 
eventually  with  the  model. 

The  model  inputs  for  the  meso-scale  laboratory  reactor  test  are  shown  in  Figure  12.  The  reactor  length 
was  set  to  0.2 1  m  or  slightly  greater  than  the  distance  to  the  first  sampling  port.  Default  values  were 
used  for  the  two  kinetic  reaction  parameters.  The  only  other  uncertain  inputs  were  the  two  parameters 
describing  dual-domain  mass  transfer.  These  parameters  were  varied  over  multiple  model  runs,  and 
model  results  were  compared  with  the  observed  data  for  RDX  at  the  first  sampling  port.  It  was 
detennined  that  model  results  close  to  those  observed  could  only  be  obtained  for  the  fraction  of  mobile 
domain  (/)  set  to  1.0.  With  any  immobile  domain,  the  computed  steady-state  RDX  concentrations 
exceeded  those  observed  by  as  much  as  double  those  observed  in  some  cases.  The  reaction  rate  in  the 
immobile  domain  is  much  less  than  for  the  mobile  domain  due  to  more  limiting  mass  transfer.  Dual¬ 
domain  mass  transfer  may  not  be  an  appropriate  mechanism  for  describing  the  flow  within  the  reactor. 
It  is  possible  that  the  flow  can  totally  bypass  some  regions  of  the  reactor  without  providing  any 
opportunity  for  mass  transfer  from  flowing  water  to  dead-end  pore  spaces. 


"Meso-scale  application  with  ID 

surface-site  reactor  m 

"Length,  m" 

0.21 

"Width,  m" 

.508 

"Height,  m" 

.305 

"Porosity" 

0.50 

"fraction  mobile" 

1. 

"Dispersivity  constant" 

0.1 

"domain  mass  trans  rate,  1/hr" 

0.1 

"RDX  mole,  diffus.,  cm2/sec" 

5.9E-6 

"packing  solids  diam,  m" 

.01 

"Fe  mass  in  reactor,  kg" 

2.0 

"Fe  half  sat  cone,  kO,  mol/L" 

1.806 

"Fe  max  reac  rate,  kmax,  1/min" 

2.824 

"RDX  initial,  mg/L" 

0.0 

"number  of  delta  x" 

10 

"delta  t,  hr" 

.002 

"simulation  duration,  hr" 

2.0 

"DX  output  intervals" 

1 

"solution  scheme" 

1 

"Time,  hr  flowin,  L/hr  RDX  in. 

mg/L  Tracer  in,  mg/L" 

0.0  120.0 

.77  1.0 

2.0  120.0 

.77  1.0 

Figure  12. Input  parameters  for  longitudinal  reactor 
model,  meso-scale  laboratory  reactor  test. 


Model-computed  RDX  concentrations  at  X  =  0.2  m  for  the  inputs  shown  in  Figure  12  are  plotted 
versus  time  in  Figure  13  along  with  observed  RDX  concentrations  at  that  location.  Again,  this  is  for/= 
1.0.  The  model  compares  quite  well  with  the  observations,  even  when  using  the  default  kinetic  reaction 
parameters.  Slight  adjustments  to  these  parameters  can  improve  agreement  with  observed  results. 
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Figure  ^.Concentrations  of  RDX  versus  time  at  X  = 

0.2  m  for  first  sampling  port  of  the  meso-scale 
laboratory  reactor. 

CONCLUSIONS:  This  model  formulation  includes  the  coupling  of  mass  transfer  of  dissolved 
constituents  from  pore  water  to  iron  solid  surfaces  and  surface  site  reduction-reaction  kinetics.  Two 
model  versions  were  developed,  a  batch  reactor  version,  and  a  ID,  longitudinal  reactor  version.  The 
latter  version  is  far  more  versatile  for  simulating  field  conditions  since  it  can  handle  time-varying 
influent  loading  and  water  flow  rate.  The  good  agreement  between  model  results  and  meso-scale 
laboratory  data  indicates  that  the  default  kinetic  reaction  parameters  may  be  applicable,  and  the  use  of 
dual-domain  mass  transfer  may  be  inappropriate  and  not  needed.  Additional  laboratory  and/or  field 
data  are  required  to  complete  model  testing,  evaluation,  and  validation. 
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